Circular and non-circular nearly horizon-skimming orbits in Kerr spacetimes 
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We have performed a detailed analysis of orbital motion in the vicinity of a nearly extremal Kerr black hole. 
For very rapidly rotating black holes — spin parameter a = J/M > 0.9524M — we have found a class of 
very strong field eccentric orbits whose orbital angular momentum L z increases with the orbit's inclination with 
respect to the equatorial plane, while keeping latus rectum and eccentricity fixed. This behavior is in contrast 
with Newtonian intuition, and is in fact opposite to the "normal" behavior of black hole orbits. Such behavior 
was noted previously for circular orbits; since it only applies to orbits very close to the black hole, they were 
named "nearly horizon-skimming orbits". Our current analysis generalizes this result, mapping out the full 
generic (inclined and eccentric) family of nearly horizon- skimming orbits. The earlier work on circular orbits 
reported that, under gravitational radiation emission, nearly horizon- skimming orbits exhibit unusual inspiral, 
tending to evolve to smaller orbit inclination, toward prograde equatorial configuration. Normal orbits, by 
contrast, always demonstrate slowly growing orbit inclination — orbits evolve toward the retrograde equatorial 
configuration. Using up-to-date Teukolsky-based fluxes, we have concluded that the earlier result was incorrect 
— all circular orbits, including nearly horizon-skimming ones, exhibit growing orbit inclination under radiative 
backreaction. Using kludge fluxes based on a Post-Newtonian expansion corrected with fits to circular and 
to equatorial Teukolsky-based fluxes, we argue that the inclination grows also for eccentric nearly horizon- 
skimming orbits. We also find that the inclination change is, in any case, very small. As such, we conclude that 
these orbits are not likely to have a clear and peculiar imprint on the gravitational waveforms expected to be 
measured by the space-based detector LISA. 
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I. INTRODUCTION 

The space-based gravitational-wave detector LISA |[l|] will 
be a unique tool to probe the nature of supermassive black 
holes (SMBHs), making it possible to map in detail their 
spacetimes. This goal is expected to be achieved by observing 
gravitational waves emitted by compact stars or black holes 
with masses fj, as 1 — 100 Af© spiraling into the SMBHs 
which reside in the center of galaxies [2] (particularly the low 
end of the galactic center black hole mass function, M « 
10 5 — 10 7 Mq). Such events are known as extreme mass ratio 
inspirals (EMRIs). Current wisdom suggests that several tens 
to perhaps of order a thousand such events could be measured 
per year by LISA Jl. 

Though the distribution of spins for observed astrophysi- 
cal black holes is not very well understood at present, very 
rapid spin is certainly plausible, as accretion tends to spin- 
up SMBHs [4]. Most models for quasi-periodic oscillations 
(QPOs) suggest this is indeed the case in all low-mass x-ray 
binaries for which data is available [5|]. On the other hand, 
continuum spectral fitting of some high-mass x-ray binaries 
indicates that modest spins (spin parameter a/M = J/M 2 ~ 
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0.6 — 0.8) are likewise plausible [6]. The continuum-fit tech- 
nique does find an extremely high spin of a/M > 0.98 for the 
galactic "microquasar" GRS1915+105 [7]. This argues for a 
wide variety of possible spins, depending on the detailed birth 
and growth history of a given black hole. 

In the mass range corresponding to black holes in galactic 
centers, measurements of the broad iron Ka emission line in 
active galactic nuclei suggest that SMBHs can be very rapidly 
rotating (see Ref. [8] for a recent review). For instance, in 
the case of MCG-6-30-15, for which highly accurate observa- 
tions are available, a has been found to be larger than 0.987M 
at 90% confidence [9]. Because gravitational waves from EM- 
RIs are expected to yield a very precise determination of the 
spins of SMBHs 111 Oil , it is interesting to investigate whether 
EMRIs around very rapidly rotating black holes may possess 
peculiar features which would be observable by LISA. Should 
such features exist, they would provide unambiguous infor- 
mation on the spin of SMBHs and thus on the mechanisms 
leading to their formation flUl . 

For extremal Kerr black holes (a — M), the existence of 
a special class of "circular" orbits was pointed out long ago 
by Wilkins 11211 . who named them "horizon-skimming" or- 
bits. ("Circular" here means that the orbits are of constant 
Boyer-Lindquist coordinate radius r.) These orbits have vary- 
ing inclination angle with respect to the equatorial plane and 
have the same coordinate radius as the horizon, r = M. De- 
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spite this seemingly hazardous location, it can be shown that 
all these r — M orbits have finite separation from one another 
and from the event horizon [13]. Their somewhat pathological 
description is due to a singularity in the Boyer-Lindquist co- 
ordinates, which collapses a finite span of the spacetime into 
r = M. 

Besides being circular and "horizon-skimming", these or- 
bits also show peculiar behavior in their relation of angular 
momentum to inclination. In Newtonian gravity, a generic or- 
bit has L z = \L\ cos i, where i is the inclination angle relative 
to the equatorial plane (going from i — for equatorial pro- 
grade orbits to i = w for equatorial retrograde orbits, passing 
through l — ir/2 for polar orbits), and L is the orbital angular 
momentum vector. As a result, dL z (r, t) /di < 0, and the an- 
gular momentum in the z-direction always decreases with in- 
creasing inclination if the orbit's radius is kept constant. This 
intuitively reasonable decrease of L z with l is seen for almost 
all black hole orbits as well. Horizon-skimming orbits, by 
contrast, exhibit exactly the opposite behavior: L z increases 
with inclination angle. 

Reference Hill asked whether the behavior dL z /di > 
could be extended to a broader class of circular orbits than just 
those at the radius r = M for the spin value a — M. It was 
found that this condition is indeed more general, and extended 
over a range of radius from the "innermost stable circular or- 
bit" to r ~ 1.8M for black holes with a > 0.9524M. Or- 
bits that show this property have been named "nearly horizon- 
skimming". The Newtonian behavior dL z (r, C)jdi < is 
recovered for all orbits at r > 1.8M lfl4ll . 

A qualitative understanding of this behavior comes from 
recalling that very close to the black hole all physical pro- 
cesses become "locked" to the hole's event horizon lfl5ll . with 
the orbital motion of point particles coupling to the horizon's 
spin. This locking dominates the "Keplerian" tendency of an 
orbit to move more quickly at smaller radii, forcing an or- 
biting particle to slow down in the innermost orbits. Lock- 
ing is particularly strong for the most-bound (equatorial) or- 
bits; the least-bound orbits (which have the largest inclination) 
do not strongly lock to the black hole's spin until they have 
very nearly reached the innermost orbit 11411 . The property 
dL z (r, i)/di > reflects the different efficiency of nearly 
horizon-skimming orbits to lock with the horizon. 

Reference [14] argued that this behavior could have ob- 
servational consequences. It is well-known that the incli- 
nation angle of an inspiraling body generally increases due 
to gravitational-wave emission Ilia 1 1711 . Since dL z /dt < 
because of the positive angular momentum carried away 
by the gravitational waves, and since "normal" orbits have 
dL z /dt < 0, one would expect di/dt > 0. However, if 
during an evolution dL z /di switches sign, then di/dt might 
switch sign as well: An inspiraling body could evolve towards 
an equatorial orbit, signalling the presence of an "almost- 
extremal" Kerr black hole. 

It should be emphasized that this argument is not rigorous. 
In particular, one needs to consider the joint evolution of or- 
bital radius and inclination angle; and, one must include the 
dependence of these two quantities on orbital energy as well 



as angular momentum 1 . As such, di/dt depends not only on 

dL z /dt and dL z /di, but also on dE/dt, dE/di, dE/dr and 
dL z /dr. 

In this sense, the argument made in Ref. [14] should be 
seen as claiming that the contribution coming from dL z /dt 
and dL z /dt are simply the dominant ones. Using the nu- 
merical code described in [ 17] to compute the fluxes dL z /dt 
and dE/dt, it was then found that a test-particle on a circu- 
lar orbit passing through the nearly horizon-skimming region 
of a Kerr black hole with a = 0.998A/ (the value at which 
a hole's spin tends to be buffered due to photon capture from 
thin disk accretion IT19I1 ) had its inclination angle decreased 
by Si w 1° — 2° H14fl in the adiabatic approximation 12011 . 
It should be noted at this point that the rate of change of in- 
clination angle, di/dt, appears as the difference of two rel- 
atively small and expensive to compute rates of change [cf. 
Eq. (3.8) of Ref. [ 17]]. As such, small relative errors in those 
rates of change can lead to large relative errors in di/dt. Fi- 
nally, in Ref. lfl4ll it was speculated that the decrease could 
be even larger for eccentric orbits satisfying the condition 
dL z /di > 0, possibly leading to an observable imprint on 
EMRI gravitational waveforms. 

The main purpose of this paper is to extend Ref. B14ll 's anal- 
ysis of nearly horizon-skimming orbits to include the effect 
of orbital eccentricity, and to thereby test the speculation that 
there may be an observable imprint on EMRI waveforms of 
nearly horizon-skimming behavior. In doing so, we have re- 
visited all the calculations of Ref. iflill using a more accurate 
Teukolsky solver which serves as the engine for the analysis 
presented in Ref. 12111 . 

We have found that the critical spin value for circular 
nearly horizon-skimming orbits, a > 0.9524M, also delin- 
eates a family of eccentric orbits for which the condition 
dL z (p, e, i)/di > holds. (More precisely, we consider vari- 
ation with respect to an angle 0- mc that is easier to work with 
in the extreme strong field, but that is easily related to t.) The 
parameters p and e are the orbit's latus rectum and eccentric- 
ity, defined precisely in Sec. [II] These generic nearly horizon- 
skimming orbits all have p < 2M, deep in the black hole's 
extreme strong field. 

We next study the evolution of these orbits under 
gravitational-wave emission in the adiabatic approximation. 
We first revisited the evolution of circular, nearly horizon- 
skimming orbits using the improved Teukolsky solver which 
was used for the analysis of Ref. 12111 . The results of this anal- 
ysis were somewhat surprising: Just as for "normal" orbits, 
we found that orbital inclination always increases during in- 
spiral, even in the nearly horizon-skimming regime. This is in 
stark contrast to the claims of Ref. [ 14]. As noted above, the 
inclination's rate of change depends on the difference of two 
expensive and difficult to compute numbers, and thus can be 
strongly impacted by small relative errors in those numbers. 



In the general case, one must also include the dependence on "Carter's 
constant" Q II 811 . the third integral of black hole orbits (described more 
carefully in Sec.UJJ. For circular orbits, Q = Q(E, L z ): knowledge of E 
and L z completely determines Q. 



3 



A primary result of this paper is thus to retract the claim of 
Ref. tl4j that an important dynamic signature of the nearly 
horizon-skimming region is a reversal in the sign of incli- 
nation angle evolution: The inclination always grows under 
gravitational radiation emission. 

We next extended this analysis to study the evolution of 
generic nearly horizon-skimming orbits. The Teukolsky code 
to which we have direct access can, at this point, only com- 
pute the radiated fluxes of energy E and angular momentum 
L z \ results for the evolution of the Carter constant Q are just 
now beginning to be understood B22I1 . and have not yet been 
incorporated into this code. We instead use "kludge" expres- 
sions for dEJdt, dL z /dt, and dQ/dt which were inspired by 
Refs. These expression are based on post-Newtonian 

flux formulas, modified in such a way that they fit strong-field 
radiation reaction results obtained from a Teukolsky integra- 
tor; see Ref. [24] for further discussion. Our analysis indicates 
that, just as in the circular limit, the result di/dt > holds 
for generic nearly horizon-skimming orbits. Furthermore, and 
contrary to the speculation of Ref. IU4I1 . we do not find a large 
amplification of di/dt as orbits are made more eccentric. 

Our conclusion is that the nearly horizon-skimming regime, 
though an interesting curiosity of strong-field orbits of nearly 
extremal black holes, will not imprint any peculiar observa- 
tional signature on EMRI waveforms. 

The remainder of this paper is organized as follows. In Sec. 
II, we review the properties of bound stable orbits in Kerr, pro- 
viding expressions for the constants of motion which we will 
use in Sec. Ill to generalize nearly horizon-skimming orbits to 
the non-circular case. In Sec. IV, we study the evolution of the 
inclination angle for circular nearly horizon-skimming orbits 
using Teukolsky-based fluxes; in Sec. V we do the same for 
non-circular orbits and using kludge fluxes. We present and 
discuss our detailed conclusions in Sec. VI. The fits and post- 
Newtonian fluxes used for the kludge fluxes are presented in 
the Appendix. Throughout the paper we have used units in 
which G = c = 1. 



H. BOUND STABLE ORBITS IN KERR SPACETIMES 

The line element of a Kerr spacetime, written in Boyer- 
Lindquist coordinates reads lEBIl 



ds 2 
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(1) 



where 



E = r 2 + a 2 cos 2 ' 



A = r — 2M r + a. (2) 



Up to initial conditions, geodesies can then be labelled by four 
constants of motion: the mass p, of the test particle, its energy 
E and angular momentum L z as measured by an observer at 



infinity and the Carter constant Q |fl8ll . The presence of these 
four conserved quantities makes the geodesic equations sepa- 
rable in Boyer-Lindquist coordinates. Introducing the Carter 
time A, defined by 
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dX 
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the geodesic equations become 
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The conserved parameters E, L z , and Q can be remapped 
to other parameters that describe the geometry of the orbit. We 
have found it useful to describe the orbit in terms of an angle 
#min — the minimum polar angle reached by the orbit — as 
well as the latus rectum p and the eccentricity e. In the weak- 
field limit, p and e correspond exactly to the latus rectum and 
eccentricity used to describe orbits in Newtonian gravity; in 
the strong field, they are essentially just a convenient remap- 
ping of the orbit's apoastron and periastron: 
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P 



' ap 



1 ' P crl i i 

1 — e 1 + e 

Finally, in much of our analysis, it is useful to refer to 

Z_ EE cos 2 6> min , (8) 

rather than to 9 min directly. 

To map (E,L Z ,Q) to (p, e, z_), use Eq. (0]l to impose 
dr/dX = at r = r ap and r = r pci i, and to impose 
d9/dX = at 9 = 9 m [ n . (Note that for a circular orbit, 
^ap = fperi = tq. In this case, one must apply the rules 
dr/dX = and d 2 r/dX 2 = at r = r .) Following this ap- 
proach, Schmidt [26] was able to derive explicit expressions 
for E, L z and Q in terms of p, e and Z-. We now briefly 
review Schmidt's results. 

Let us first introduce the dimensionless quantities 
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Figure 1: Le/i panel: Inclination angles 9 1Ilc for which bound stable orbits exist for a black hole with spin a — 0.998 M. The allowed range 
for 6*i nc goes from 6i nc = to the curve corresponding to the eccentricity under consideration, 6i nc = 6™*- Right panel: Same as left but for 
an extremal black hole, a = M. Note that in this case never reaches zero. 
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and the functions 
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Let us further define the set of functions 
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The energy of the particle can then be written 
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The parameter D takes the values ±1. The angular momen- 
tum is a solution of the system 

hE 2 -2g 1 EL z ~h 1 L 2 z -d 1 = 0, (23) 
f 2 E 2 ~2g 2 EL z -h 2 L 2 z -d 2 =0. (24) 

By eliminating the L 2 terms in these equations, one finds the 
solution 



pE 2 - K 
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(25) 



for the angular momentum. Using d9/d\ = at 8 = 8„ 
the Carter constant can be written 
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Additional constraints on p, e, z_ are needed for the orbits 
to be stable. Inspection of Eq. shows that an eccentric orbit 
is stable only if 



dV r 
dr 



(Vpcri) > . 



(27) 



It is marginally stable if dV r /dr = at r — r per i. Similarly, 
the stability condition for circular orbits is 



d 2 V r 
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marginally stable orbits are set by d 2 V r /dr 2 = at r = Tq. 
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Finally, we note that one can massage the above solutions 
for the conserved orbital quantities of bound stable orbits to 
rewrite the solution for L z as 



L, 



giE 
fa 



D 
fa 



(fiE 2 - fa)fa 



(29) 



From this solution, we see that it is quite natural to refer to 
orbits with D = 1 as prograde and to orbits with D = — 1 
as retrograde. Note also that Eq. (l29l is a more useful form 
than the corresponding expression, Eq. (A4), of Ref. 112 ill . In 
that expression, the factor 1/fa has been squared and moved 
inside the square root. This obscures the fact that hi changes 
sign for very strong field orbits. Differences between Eq. (|29l 
and Eq. (A4) of [21] are apparent for a > 0.835, although 
only for orbits close to the separatrix (i.e., the surface in the 
space of parameters (p, e, i) where marginally stable bound 
orbits lie). 



III. NON-CIRCULAR NEARLY HORIZON-SKIMMING 
ORBITS 

With explicit expressions for E, L z and Q as functions of 
p, e and z~, we now examine how to generalize the condi- 
tion dL z (r, l) /di > 0, which defined circular nearly horizon- 
skimming orbits in Ref. 111411 . to encompass the non-circular 
case. We recall that the inclination angle u is defined as H 1 411 



L, 



(30) 



Such a definition is not always easy to handle in the case of 
eccentric orbits. In addition, i does not have an obvious phys- 
ical interpretation (even in the circular limit), but rather was 
introduced essentially to generalize (at least formally) the def- 
inition of inclination for Schwarzschild black hole orbits. In 
that case, one has Q = L 2 . + L 2 and therefore L z = \L\ cos i. 

A more useful definition for the inclination angle in Kerr 
was introduced in Ref. Q2lll: 



D8„ 



(31) 



angle i, particularly the "kludge" computation of fluxes de- 
scribed in Sec.[V] (This is because the kludge fluxes are based 
on an extension of post-Newtonian formulas to the strong- 
field regime, and these formulas use l for inclination angle.) 
Accordingly, we often switch back and forth between these 
two notions of inclination, and in fact present our final results 
for inclination evolution using both du / dt and d9i nc / dt. 

Before mapping out the region corresponding to nearly 
horizon-skimming orbits, it is useful to examine stable or- 
bits more generally in the strong field of rapidly rotating black 
holes. We first fix a value for a, and then discretize the space 
of parameters (p, e, 0- mc ). We next identify the points in this 
space corresponding to bound stable geodesic orbits. Suffi- 
ciently close to the horizon, the bound stable orbits with spec- 
ified values of p and e have an inclination angle 9 inc ranging 
from (equatorial orbit) to a maximum value 6^ x . For given 
p and e, flFJO* defines the separatrix between stable and unsta- 
ble orbits. 

Example separatrices are shown in Fig.Q]for a — 0.998M 
and a = M. This figure shows the behavior of 9™ x as a 
function of the latus rectum for the different values of the ec- 
centricity indicated by the labels. Note that for a = 0.998A/ 
the angle 0f™ x eventually goes to zero. This is the general 
behavior for a < M. On the other hand, for an extremal black 
hole, a = M, 9"^ x never goes to zero. The orbits which re- 
side at r = M (the circular limit) are the "horizon-skimming 
orbits" identified by Wilkins B12ll ; the a = M separatrix has 
a similar shape even for eccentric orbits. As expected, we 
find that for given latus rectum and eccentricity the orbit with 
#inc = is the one with the lowest energy E (and hence is the 
most-bound orbit), whereas the orbit with 9 mc = 9\™ x has the 
highest E (and is least bound). 

Having mapped out stable orbits in (p, e, 9i nc ) space, we 
then computed the partial derivative dL z (p, e, 9i nc ) / d9 lrLC and 
identified the following three overlapping regions: 

• Region A: The portion of the (p, e) plane for which 
dL z {p,e,9 inc )/d9 inc > OforO < 8 inc < 0£» This 
region is illustrated in Fig.|2]as the area under the heavy 
solid line and to the left of the dot-dashed line (green in 
the color version). 



where 6* m i n is the minimum reached by 8 during the orbital 
motion. This angle is trivially related to z_ (z- — sin 2 9i nc ) 
and ranges from to ir/2 for prograde orbits and from ir/2 to 
7r for retrograde orbits. It is a simple numerical calculation to 
convert between i and 9 mc ; doing so shows that the differences 
between i and 8 inc are very small, with the two coinciding for 
a = 0, and with a difference that is less than 2.6° for a = M 
and circular orbits with r = M. 

Bearing all this in mind, the condition we have adopted to 
generalize nearly horizon-skimming orbits is 



dL z (p, e, 9 h 
d9; nc 



> 



(32) 



We have found that certain parts of this calculation, particular 
the analysis of strong-field geodesic orbits, are best done us- 
ing the angle 6*; nc ; other parts are more simply done using the 



• Region B: The portion of the 
for which (L z ) most bound (p, e) is 
{Lz ) least bound (p, e). In other words, 



L z (p,e,0) < L z (p, 



nmax 
inc 



(p, e) plane 
smaller than 



(33) 



in Region B. Note that Region B contains Region A. It 
is illustrated in Fig. [2] as the area under the heavy solid 
line and to the left of the dotted line (red in the color 
version). 

• Region C: The portion of the (p, e) plane for which 
dL z (p,e,9i nc )/d9 inc > for at least one angle 9{ nc 
between and 6*|"^ x . Region C contains Region B, and 
is illustrated in Fig. [2] as the area under the heavy solid 
line and to the left of the dashed line (blue in the color 
version). 
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Figure 2: Left panel: Non-circular nearly horizon-skimming orbits for a = 0.998Af. The heavy solid line indicates the separatrix between 
stable and unstable orbits for equatorial orbits (t = #i Ilc = 0). All orbits above and to the left of this line are unstable. The dot-dashed line 
(green in the color version) bounds the region of the (p, e)-plane where dL z /d6 luc > for all allowed inclination angles ("Region A"). All 
orbits between this line and the separatrix belong to Region A. The dotted line (red in the color version) bounds the region (L z ) most bound < 
(-t'z)least bound ("Region B"). Note that B includes A. The dashed line (blue in the color version) bounds the region where dL z /d8i n c > for 
at least one inclination angle ("Region C"); note that C includes B. All three of these regions are candidate generalizations of the notion of 
nearly horizon- skimming orbits. Right panel: Same as the left panel, but for the extreme spin case, a — M. In this case the separatrix between 
stable and unstable equatorial orbits is given by the line p/M = 1 + e. 



Orbits in any of these three regions give possible general- 
izations of the nearly horizon-skimming circular orbits pre- 
sented in Ref. 11411 . Notice, as illustrated in Fig. [2] that the 
size of these regions depends rather strongly on the spin of 
the black hole. All three regions disappear altogether for 
a < 0.9524A/ (in agreement with their sizes grow with 
a, reaching maximal extent for a = M . These regions never 
extend beyond p ~ 2 M. 

As we shall see, the difference between these three regions 
is not terribly important for assessing whether there is a strong 
signature of the nearly horizon-skimming regime on the inspi- 
ral dynamics. As such, it is perhaps most useful to use Region 
C as our definition, since it is the most inclusive. 



IV. EVOLUTION OF inc : CIRCULAR ORBITS 

To ascertain whether nearly horizon-skimming orbits can 
affect an EMRJ in such a way as to leave a clear imprint in the 
gravitational-wave signal, we have studied the time evolution 
of the inclination angle 6*j llc . To this purpose we have used the 
so-called adiabatic approximation H20fl . in which the infalling 
body moves along a geodesic with slowly changing parame- 
ters. The evolution of the orbital parameters is computed us- 
ing the time-averaged fluxes dE/dt, dL z /dt and dQ/dt due 
to gravitational-wave emission ("radiation reaction"). As dis- 
cussed in Sec. mi E, L z and Q can be expressed in terms of p, 



e, and f? mc . Given rates of change of E, L z and Q, it is then 
straightforward |23(] to calculate dp/dt, de/dt, and d9 mc /dt 
(or du/dt). 

We should note that although perfectly well-behaved for all 
bound stable geodesies, the adiabatic approximation breaks 
down in a small region of the orbital parameters space very 
close to the separatrix, where the transition from an inspiral to 
a plunging orbit takes place i27ll . However, since this region 
is expected to be very small 2 and its impact on LISA wave- 
forms rather hard to detect 12711 . we expect our results to be 
at least qualitatively correct also in this region of the space of 
parameters. 

Accurate calculation of dE/dt and dL z /dt in the adiabatic 
approximation involves solving the Teukolsky and Sasaki- 
Nakamura equations I281I29I1 . For generic orbits this has been 
done for the first time in Ref. l2~fll . The calculation of dQ / dt 
for generic orbits is more involved. A formula for dQ /dt has 
been recently derived [22], but has not yet been implemented 
(at least in a code to which we have access). 

On the other hand, it is well-known that a circular orbit will 
remain circular under radiation reaction lf30L l3ll l32h . This 
constraint means that Teukolsky-based fluxes for E and L z 



2 Its width in p/M is expected to be of the order of Ap/M ~ (/1/M) 2 / 5 , 
where fi is the mass of the infalling body 2711. 
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are sufficient to compute dQ/dt. Considering this limit, the 
rate of change dQ jdt can be expressed in terms of dE / dt and 

dL I dt as 



dQ 
dt 



N^p, i) (dE 

■jv b (p,0 V * 



at 5 (p, 



dt 



(34) 



where 



JVi(p,t) = s(p, Op 4 + o 2 s( p , Op 2 

-2aM(L z (p,L)-aE(p,i))p, (35) 
iV 4 (p, i) = (2 M p - p 2 ) (p, t) - 2 M a E(p, l) p , (36) 



N 5 {p,i) = (2Mp-p 2 -a 2 )/2 



(37) 



(These quantities are for a circular orbit of radius p.) Using 
this, it is simple to compute dd- mc / dt (or djjdt). 

This procedure was followed in Ref. IU4I1 . using the code 
presented in Ref. lfl7ll . to determine the evolution of i\ this 
analysis indicated that di/dt < for circular nearly horizon- 
skimming orbits. As a first step to our more general analy- 
sis, we have repeated this calculation but using the improved 
Sasaki-Nakamura-Teukolsky code presented in Ref. 02111 : we 
focused on the case a — 0.998M. 

Rather to our surprise, we discovered that the fluxes dE/dt 
and dL z /dt computed with this more accurate code indicate 
that di/dt > (and dO inc /dt > 0) for all circular nearly 
horizon-skimming orbits — in stark contrast with what was 
found in Ref. IU4I1 . As mentioned in the introduction, the rate 
of change of inclination angle appears as the difference of two 
quantities. These quantities nearly cancel (and indeed cancel 
exactly in the limit a = 0); as such, small relative errors in 
their values can lead to large relative error in the inferred in- 
clination evolution. Values for dE/dt, dL z /dt, di/dt, and 
ddinc/dt computed using the present code are shown in Ta- 
ble|I]in the columns with the header "Teukolsky". 



V. EVOLUTION OF inc : NON-CIRCULAR ORBITS 

The corrected behavior of circular nearly horizon- 
skimming orbits has naturally led us to investigate the evo- 
lution of non-circular nearly horizon-skimming orbits. Since 
our code cannot be used to compute dQ / dt, we have resorted 
to a "kludge" approach, based on those described in Refs. 
Ir23l l24ll . In particular, we mostly follow the procedure de- 
veloped by Gair & Glampedakis B2411 . though (as described 
below) importantly modified. 

The basic idea of the "kludge" procedure is to use the func- 
tional form of 2PN fluxes E, L z and Q, but to correct the 
circular part of these fluxes using fits to circular Teukolsky 
data. As developed in Ref. Il24ll . the fluxes are written 



dt 



(l_ e 2)3/2 (1 _ e .) 



GG 



.2\-3/2 



dLz 

dt 



(p, e, i) 

2PN 



dLz 
dt 



(p,0,i) + 

2PN 



dL, 
dt 



(p, 

fit circ 



(39) 



(^) G =(l-e 2 ) 3/2 vW^)x 



fdQ/dt\ 



fdQ/dt\ 



(1 _ e2) - 3/2 , r-_p I {p> e , o - I ^ ) (p, 0, i) 

\ VV / 2PN V VV / 2PN 



f dQ/dt \ 

\ VQ J fit circ 



(40) 



The post-Newtonian fluxes (dE /rft)2PN> {dL z /dt)2PN and 
(<2Q/di)2PN are given in the Appendix [particularly Eqs. 
dAll. SK5\. and (TOlll. 

Since for circular orbits the fluxes dE/dt, dL z /dt and 
dQ/dt are related through Eq. d34b , only two fits to circu- 
lar Teukolsky data are needed. One possible choice is to fit 
dL z /dt and di/dt, and then use the circularity constraint to 
obtain 3 1241 



\ VV /fit circ 

dL^ 
dt 



2 tan i 



' • 2 

fit circ Sm L 



fit circ 



(41) 



dE\ 

-jr (P, 

al / fit circ 



Njjp, i) 

'jvi(p,0 



al I fit circ 



N 5 (p, i) / nf n v dQ dt \ 

— — -VQ(P,0,i) —f=- (p,t) ■ (42) 

Nl(p,i) V VQ /fit circ 



As stressed in Ref. [24], one does not expect these fluxes to 
work well in the strong field, both because the post-Newtonian 
approximation breaks down close to the black hole, and be- 
cause the circular Teukolsky data used for the fits in Ref. ll24ll 
was computed for 3M < p < 30M. As a first attempt to 
improve their behavior in the nearly horizon-skimming re- 
gion, we have made fits using circular Teukolsky data for 
orbits with M < p < 2M. In particular, for a black hole 
with a = 0.998M, we computed the circular Teukolsky-based 
fluxes dL z /dt and di/dt listed in Table U (columns 8 and 10). 
These results were fit (with error < 0.2%); see Eqs. ( 1A.41 i and 
(1A.61 > in the Appendix. 



dE\ 
~dt) 



GG 



.2^3/2 



dE 



(1 



<§) (!>■.<•,>■) 

/ 2PN 



dt 



(p,0,0 



2PN 



dE 
~dt 



(p, 



lit c 



(38) 



3 This choice might seem more involved than fitting directly dL z /dt and 
dQ/dt, but, as noted by Gair & Glampedakis, ensures more sensible re- 
sults for the evolution of the inclination angle. This generates more physi- 
cally realistic inspirals 12411 . 
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Despite using strong-field Teukolsky fluxes for our fit, we 
found fairly poor behavior of these rates of change, particu- 
larly as a function of eccentricity. To compensate for this, we 
introduced a kludge-type fit to correct the equatorial part of 
the flux, in addition to the circular part. We fit, as a function 
of p and e, Teukolsky -based fluxes for dE/dt and dL z /dt for 
orbits in the equatorial plane, and then introduce the following 
kludge fluxes for E and L z : 



dE 

--fee,,) 



dL 



dE\ 

-jr fe e, l) 

) GG 
dE \ , 

-jr (p,e,0)H 

/ GG 

lit) (p ' e,i) 

UL J GG 

^) fee,0) 

dt / GG 



^) fee) (43) 

dt /fiteq 



^) fee). (44) 

d< / fiteq 



[Note that Eq. ( f40b for c?Q /dt is not modified by this proce- 
dure since dQ/dt = for equatorial orbits.] Using equatorial 
non-circular Teukolsky data provided by Drasco |2 ll l33ll for 
a = 0.998 and M < p < 2M (the t = "Teukolsky" data in 
Tables HI1 HITl and ITVli, we found fits (with error < 1.5%); see 
Eqs. dA.9t and dA.lQI ). Note that the fits for equatorial fluxes 
are significantly less accurate than the fits for circular fluxes. 
This appears to be due to the fact that, close to the black hole, 
many harmonics are needed in order for the Teukolsky-based 
fluxes to converge, especially for eccentric orbits (cf. Figs. 2 
and 3 of Ref. 12111 . noting the number of radial harmonics that 
have significant contribution to the flux). Truncation of these 
sums is likely a source of some error in the fluxes themselves, 
making it difficult to make a fit of as high quality as we could 
in the circular case. 

These fits were then finally used in Eqs. (l43l and (l44b to 
calculate the kludge fluxes dE/dt and dL z /dt for generic or- 
bits. This kludge reproduces to high accuracy our fits to the 
Teukolsky-based fluxes for circular orbits (e = 0) or equato- 
rial orbits (t = 0). Some residual error remains because the 
i = limit of the circular fits do not precisely equal the e = 
limit of the equatorial fits. 

Table U compares our kludge to Teukolsky-based fluxes for 
circular orbits; the two methods agree to several digits. Tables 
HTl HHI and [IV] compare our kludge to the generic Teukolsky- 
based fluxes for dE/dt and dL z /dt provided by Drasco 
112 U 13311 . In all cases, the kludge fluxes dE/dt and dL z /dt 
have the correct qualitative behavior, being negative for all 
the orbital parameters under consideration (a = 0.998A/, 
1 < p/M < 2, < e < 0.5 and 0° < i < 41°). The 
relative difference between the kludge and Teukolsky fluxes 
is always less than 25% for e = and e = 0.1 (even for 
orbits very close to separatrix). The accuracy remains good 
at larger eccentricity, though it degrades somewhat as orbits 
come close to the separatrix. 

Tables HI [TT1 iHll andllVI also present the kludge values of the 
fluxes di/dt and d0 mc /dt as computed using Eqs. d43l and 
d44l for dE/dt and dL z /dt, plus Eq. (gDJ for dQ/dt. Though 
certainly not the last word on inclination evolution (pending 



rigorous computation of dQ / dt), these rates of change proba- 
bly represent a better approximation than results published to 
date in the literature. (Indeed, prior work has often used the 
crude approximation di/dt = lf2lll to estimate dQ/dt given 
dE/dt and dL z /dt.) 

Most significantly, we find that (c?i/dt)kiudgc > and 
(dOinc/dt) kludge > for all of the orbital parameters we con- 
sider. In other words, we find that di/dt and d6- mc /dt do not 
ever change sign. 

Finally, in Table [Vl we compute the changes in 6 lnc and i 
for the inspiral with mass ratio /.i/M = 10~ 6 . In all cases, 
we start at p/M = 1.9. The small body then inspirals through 
the nearly horizon-skimming region until it reaches the sep- 
aratrix; at this point, the small body will fall into the large 
black hole on a dynamical timescale ~ M, so we terminate 
the calculation. The evolution of circular orbits is computed 
using our fits to the circular- Teukolsky fluxes of E and L z ; 
for eccentric orbits we use the kludge fluxes d40l ), d43l and 
(f44b . As this exercise demonstrates, the change in inclination 
during inspiral is never larger than a few degrees. Not only is 
there no unique sign change in the nearly horizon-skimming 
region, but the magnitude of the inclination change remains 
puny. This leaves little room for the possibility that this class 
of orbits may have a clear observational imprint on the EMRI- 
waveforms to be detected by LISA. 



VI. CONCLUSIONS 

We have performed a detailed analysis of the orbital mo- 
tion near the horizon of near-extremal Kerr black holes. We 
have demonstrated the existence of a class of orbits, which we 
have named "non-circular nearly horizon-skimming orbits", 
for which the angular momentum L z increases with the or- 
bit's inclination, while keeping latus rectum and eccentricity 
fixed. This behavior, in stark contrast to that of Newtonian 
orbits, generalizes earlier results for circular orbits lfl4ll . 

Furthermore, to assess whether this class of orbits can pro- 
duce a unique imprint on EMRI waveforms (an important 
source for future LISA observations), we have studied, in 
the adiabatic approximation, the radiative evolution of incli- 
nation angle for a small body orbiting in the nearly horizon- 
skimming region. For circular orbits, we have re-examined 
the analysis of Ref. Il4ll using an improved code for comput- 
ing Teukolsky-based fluxes of the energy and angular momen- 
tum. Significantly correcting Ref. tl4J]'s results, we found no 
decrease in the orbit's inclination angle. Inclination always 
increases during inspiral. 

We next carried out such an analysis for eccentric nearly 
horizon-skimming orbits. In this case, we used "kludge" 
fluxes to evolve the constants of motion E, L z and Q 112411 . 
We find that these fluxes are fairly accurate when compared 
with the available Teukolsky-based fluxes, indicating that they 
should provide at least qualitatively correct information re- 
garding inclination evolution. As for circular orbits, we find 
that the orbit's inclination never decreases. For both circular 
and non-circular configurations, we find that the magnitude of 
the inclination change is quite paltry — only a few degrees at 
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most. 

Quite generically, therefore, we found that the inclination 
angle of both circular and eccentric nearly horizon-skimming 
orbits never decreases during the inspiral. Revising the results 
obtained in Ref. 11411 . we thus conclude that such orbits are 
not likely to yield a peculiar, unique imprint on the EMRI- 
waveforms detectable by LISA. 
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Appendix 

In this Appendix we report the expressions for the post- 
Newtonian fluxes and the fits to the Teukolsky data necessary 
to compute the kludge fluxes introduced in SecfV] In partic- 
ular the 2PN fluxes are given by ll24Tl 
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where fi is the mass of the infalling body and where the various e-dependent coefficients are 
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The fits to the circular- Teukolsky data of Table U are instead given by 
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where 

£(s) = ai + &iar 1/2 + cjjs- 1 , i=l,...,8, 
and the numerical coefficients are given by 

4 = -278.9387 , 6^ = 84.1414, a| = 8.6679 , b\ = -9.2401 , A = -18.3362 , 5 = 24.9034 
and by the following table 
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Note that the functional form of these fits was obtained from Eqs. (57) and (58) of Ref. M24I1 by setting a (i.e., q in their 
notation) to 1. Finally, we give expressions for the fits to the equatorial Teukolsky data of tables [II] [III] and [IV] (data with i = 0, 
columns with header "Teukolsky"): 
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2 



Table I: Teukolsky-based fluxes and kludge fluxes [computed using Eqs. J40b . J43 b and J44b l for circular orbits about a hole with a = 0.998M; 
H represents the mass of the infalling body. The Teukolsky-based fluxes have an accuracy of 10" 6 . 
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Table II: As in TableUbut for non-circular orbits; the Teukolsky-based fluxes for E and L z have an accuracy of 10 -3 . Note that our code, as 
all the Teukolsky-based code that we are aware of, presently does not have the capability to compute inclination angle evolution for generic 
orbits. 
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0.2 


16 


17.1064 


-8.171x10" 


2 


-6.995x10" 


2 


-2.398x10" 


1 


-2.096x10" 


1 


7.441x10" 


2 


AT /I w 1 A - 

O.Uz4X ID 


2 


1.7 


0.2 


24 


25.3085 


-9.948x10" 


2 


-7.443x10" 


2 


-2.806x10" 


1 


-2.178x10" 


1 


1.009x10" 


1 


O OAAw 1 A - 

O.ZVUX 1U 


2 


1.7 


0.2 


32 


33.2037 


-1.257x10" 


-1 


-8.558x10" 


2 


-3.371x10" 


1 


-2.366x10" 


1 


1.175x10" 


1 


A OA/C \, 1 A - 

y.oUox 10 


2 


1.8 


0.2 








-5.965x10" 


2 


-5.962x10" 


2 


-1.927x10" 


1 


-1.926x10" 


1 







U 




1.8 


0.2 


8 


8.5789 


-6.211x10" 


-2 


-5.997x10" 


2 


-1.990x10" 


1 


-1.930x10" 


1 


2.919x10" 


2 


T 1AA\/ 1 A - 
Z.JUUX 1U 


2 


1.8 


0.2 


16 


17.0211 


-6.953x10" 


2 


-6.147x10" 


2 


-2.175x10" 


1 


-1.954x10" 


1 


5.504x10" 


2 


/I OOAv/ 1 A - 

4.J5UX 1U 


2 


1.8 


0.2 


24 


25.2283 


-8.216x10" 


-2 


-6.502x10" 


2 


-2.474x10" 


1 


-2.016x10" 


1 


7.515x10" 


2 


n^o w 1 a - 
O.UOOX 10 


2 


1.8 


0.2 


32 


33.1656 


-1.009x10" 


1 


-7.410x10" 


2 


-2.890x10" 


1 


-2.190x10" 


1 


8.839x10" 


2 


H TCO\/ 1 A - 
/.ZJO X 10 


2 


1.9 


0.2 








-5.218x10" 


-2 


-5.210x10" 


2 


-1.786x10" 


1 


-1.783x10" 


1 







U 




1.9 


0.2 


8 


8.5312 


-5.394x10" 


2 


-5.244x10" 


2 


-1.833x10" 


1 


-1.787x10" 


1 


2.197x10" 


2 


1 HC\A \s 1 A - 

1. /U4X 10 


2 


1.9 


0.2 


16 


16.9441 


-5.928x10" 


-2 


— 5.373x 10 


2 


-1.970x10" 


1 


-1.807x10" 


1 


4.156x10" 


2 


1 \ , I A - 

j.2j4X 10 


2 


1.9 


0.2 


24 


25.1518 


-6.843x10" 


2 


-5.669x10" 


2 


-2.192x10" 


1 


-1.858x10" 


1 


5.706x10" 


2 


J J J X 1U 


2 


1.9 


0.2 


32 


33.1207 


-8.213x10" 


-2 


— 6.277x 10" 


2 


-2.502x10" 


1 


-1.966x10" 


1 


6.767x10" 


2 


C A>JC w 1 A~ 


2 


2.0 


0.2 








-4.528x10" 


2 


— 4.530x 10" 


2 


-1.637x10" 




-1.638x10" 


;[ 







A 

u 




2.0 


0.2 


8 


8.4891 


-4.657x10" 




— 4.557x 10" 


2 


-1.671x10" 




-1.641x10" 




1.679x10" 




1 .Zoj X 1U 


2 


2.0 


0.2 


16 


16.8749 


-5.049x10" 




—4.664x10" 


2 


-1.774x10" 




-1.657x10" 




3.184x10" 




Z.4j / X 1U 


2 


2.0 


0.2 


24 


25.0802 


-5.725x10" 




—4.904x10" 


2 


-1.941x10" 




-1.696x10" 




4.391x10" 




a Aw i a- 
0.44U X 1U 


2 


2.0 


0.2 


32 


33.0730 


-6.743x10" 




— 5.427x 10" 


2 


-2.175x10" 




-1.793x10" 




5.243x10" 




4. 154 X 1U 


2 


1.5 


0.3 








-8.481x10" 


2 


-8.478x10" 


~2~ 


-2.094x10" 


1 


-2.094x10" 


1 







U 




1.5 


0.3 


8 


8.7037 


-1.006x10" 


1 


-7.824x10" 


2 


-2.442x10" 


1 


-1.934x10" 


1 


1.484x10" 


1 


1 Q A 1 \/ 1 A - 

1.3U1 X 10 


1 


1.5 


0.3 


16 


17.2003 


-1.469x10" 


-1 


-7.811x10" 


2 


-3.435x10" 


1 


-1.864x10" 


1 


2.766x10" 


1 


T /I /I A w 1 A - 

Z.44UX 1U 


1 


1.6 


0.3 








-8.144x10" 


2 


-8.123x10" 


2 


-2.183x10" 


1 


-2.178x10" 


1 







U 




1.6 


0.3 


8 


8.6498 


-9.182x10" 


-2 


-7.807x10" 


2 


-2.426x10" 


1 


-2.095x10" 


1 


1.028x10" 


1 


o ni o w 1 a - 
o.91o X ID 


2 


1.6 


0.3 


16 


17.1246 


-1.223x10" 


1 


-8.089x10" 


2 


-3.122x10" 


1 


-2.144x10" 


1 


1.928x10" 


1 


1 £C3 w 1 A - 

l.ooJ X 10 


1 


1.6 


0.3 


24 


25.3046 


-1.716x10" 


-1 


-8.666x10" 


2 


-4.197x10" 


1 


-2.229x10" 


1 


2.607x10" 


1 


T OAC w 1 A - 

X 10 


1 


1.7 


0.3 








-7.362x10" 




-7.314x10" 


2 


-2.104x10" 




-2.095x10" 









U 




1.7 


0.3 


8 


8.5953 


-8.060x10" 




-7.224x10" 


2 


-2.277x10" 




-2.065x10" 




7.240x10" 




A OTyl w 1 A - 

O.ZZ4X ID 


2 


1.7 


0.3 


16 


17.0415 


-1.013x10" 




-7.369x10" 


2 


-2.774x10" 




-2.084x10" 




1.365x10" 




1 1 OA>/ 1 A - 

1. loUX 10 


1 


1.7 


0.3 


24 


25.2339 


-1.349x10" 




-7.800x10" 


2 


-3.547x10" 




-2.153x10" 




1.861x10" 




1 <COO w 1 A - 

l.OZZX ID 


1 


1.8 


0.3 








-6.488x10" 




-6.484x10" 


2 


-1.973x10" 




-1.972x10" 









U 




1.8 


0.3 


8 


8.5454 


-6.970x10" 




-6.480x10" 


2 


-2.099x10" 




-1.966x10" 




5.206x10" 




A A IKL w 1 A - 
4.4J0X 1U 


2 


1.8 


0.3 


16 


16.9628 


-8.402x10" 




-6.671x10" 


2 


-2.461x10" 




-1.998x10" 




9.857x10" 




O /I /I C w 1 A - 

o.44j X 10 


2 


1.8 


0.3 


24 


25.1601 


-1.075x10" 




-7.030x10" 


2 


-3.026x10" 




-2.056x10" 




1.353x10" 




1 1 /^n w 1 a - 
1 . 1 o9 X 1 U 


1 


1.8 


0.3 


32 


33.1047 


-1.404x10" 




-8.153x10" 


2 


-3.762x10" 




-2.255x10" 




1.600x10" 




1 QAyl 1 A - 

1.3V4X 10 


1 


1.9 


0.3 








-5.669x10" 




-5.690x10" 


2 


-1.829x10" 




-1.832x10" 









U 




1.9 


0.3 


8 


8.5010 


-6.010x10" 




-5.683x10" 


2 


-1.922x10" 




-1.824x10" 




3.823x10" 




Q T1A \/ 1 A - 

J.ZZ9X 10 


2 


1.9 


0.3 


16 


16.8911 


-7.025x10" 




— 5.8l8x 10 


2 


-2.189x10" 




-1.844x10" 




7.263x10" 




1 w 1 a - 
O.lODX 10 


2 


1.9 


0.3 


24 


25.0887 


-8.701x10" 




-6.054x10" 


2 


-2.609x10" 




-1.874x10" 




1.003x10" 




O.J / 7AiU 


2 


1.9 


0.3 


32 


33.0624 


-1.106x10" 




-6.912x10" 


2 


-3.157x10" 




-2.034x10" 




1.195x10" 




1.032x10" 


1 


2.0 


0.3 








-4.953x10" 


-2 


-4.946x10" 


2 


-1.683x10" 




-1.683x10" 


I 












2.0 


0.3 


8 


8.4616 


-5.199x10" 


2 


-4.970x10" 


2 


-1.753x10" 


1 


-1.685x10" 


1 


2.862x10" 


2 


2.395x10" 


2 


2.0 


0.3 


16 


16.8262 


-5.932x10" 


-2 


-5.079x10" 


2 


-1.954x10" 


1 


-1.699x10" 


1 


5.452x10" 


2 


4.585x10" 


2 


2.0 


0.3 


24 


25.0215 


-7.150x10" 


2 


-5.328x10" 


2 


-2.269x10" 


1 


-1.737x10" 


1 


7.564x10" 


2 


6.411x10" 


2 


2.0 


0.3 


32 


33.0172 


-8.878x10" 


-2 


-6.003x10" 


2 


-2.682x10" 


1 


-1.864x10" 


1 


9.077x10" 


2 


7.771x10" 


2 



: As in Table ITU but for additional values of eccentricity e; the Teukolsky-based fluxes for E and L z have an accuracy of 10 



V 
M 


e 


$inc 


L 


dE M 2 
dt [i z 




dE M 2 
dt i_i 2 








d£ x 7^" 




dt M 

dt x 




^inc v A/ 
dt X 






(dee ) 


(dee ) 


(kludge) 




(Teukolsky ) 






(Teukolsky] 




fkludcre^ 




(kludge) 




1.6 


0.4 








-7.766x10" 


2 


-7.772x10" 


2 


-1.918x10" 


1 


-1.919x10" 


1 







U 




1.6 


0.4 


8 


8.5863 


-9.433x10" 


2 


-7.645x10" 


2 


-2.297x10" 


1 


-1.881x10" 


1 


1.528x10" 


1 


1.3 /UX 10 


i 


1.6 


0.4 


16 


17.0151 


-1.432x10" 


1 


-7.651x10" 


2 


-3.382x10" 


1 


-1.837x10" 


1 


2.873x10" 


1 


T CO/1 \ , 1 f\ — 

z.jo4X 10 


i 


1.7 


0.4 








-7.882x10" 


2 


-7.953x10" 


2 


-2.097x10" 


1 


-2.115x10" 


1 







U 




1.7 


0.4 


8 


8.5426 


-9.002x10" 


2 


-7.408x10" 


2 


-2.367x10" 


1 


-1.978x10" 


1 


1.087x10" 


1 


y.ojox id 


2 


1.7 


0.4 


16 


16.9502 


-1.229x10" 


1 


-7.682x10" 


2 


-3.143x10" 


1 


-2.025x10" 


1 


2.054x10" 


1 


l.oJUX 1U 


1 


1.7 


0.4 


24 


25.1282 


-1.760x10" 


1 


-8.090x10" 


2 


-4.336x10" 


1 


-2.075x10" 


1 


2.809x10" 


1 


T A w 1 A - 

2.M4X 10 


1 


1.8 


0.4 








-7.107x10" 


2 


-7.007x10" 


2 


-2.013x10" 


1 


-1.988x10" 


1 







U 




1.8 


0.4 


8 


8.4989 


-7.877x10" 


2 


-7.001x10" 


2 


-2.209x10" 


1 


-1.981x10" 


1 


7.788x10" 


2 


o.o /9x 1U 


2 


1.8 


0.4 


16 


16.8817 


-1.015x10" 


1 


-7.009x10" 


2 


-2.774x10" 


1 


-1.965x10" 


1 


1.478x10" 


1 


1.3U9X 10 


1 


1.8 


0.4 


24 


25.0646 


-1.383x10" 


1 


-7.314x10" 


2 


-3.646x10" 


1 


-2.003x10" 


1 


2.036x10" 


1 


l.olUX 1U 


1 


1.8 


0.4 


32 


33.0184 


-1.887x10" 


1 


-9.193x10" 


2 


-4.755x10" 


1 


-2.319x10" 


1 


2.414x10" 


1 


Z.l JOX 10 


1 


1.9 


0.4 








-6.187x10" 


2 


-6.267x10" 


2 


-1.861x10" 


1 


-1.881x10" 


1 







U 




1.9 


0.4 


8 


8.4591 


-6.728x10" 


2 


-6.216x10" 


2 


-2.006x10" 


1 


-1.861x10" 


1 


5.666x10" 


2 


4.9oUX 1U 


2 


1.9 


0.4 


16 


16.8173 


-8.328x10" 


2 


— 6.222x 10 


2 


-2.424x10" 


1 


-1.844x10" 


1 


1.079x10" 


1 


y.jUox 1(J 


2 


1.9 


0.4 


24 


25.0006 


-1.094x10" 


1 


-6.486x10" 


2 


-3.071x10" 


1 


-1.878x10" 


1 


1.495x10" 


1 


1 in v 1 n~ 


1 


1.9 


0.4 


32 


32.9804 


-1.452x10" 


1 


-7.884x 10" 


2 


-3.896x10" 


1 


-2.158x10" 


1 


1.787x10" 


1 


1 coo \/ 1 r\ - 


1 


2.0 


0.4 








-5.483x10" 


2 


-5.457x 10" 


2 


-1.735x10" 




-1.729x10" 


;[ 







U 




2.0 


0.4 


8 


8.4235 


-5.871x10" 




—5.445x10" 


2 


-1.844x10" 




-1.720x10" 




4.222x10" 




j.OoO X LXJ 


2 


2.0 


0.4 


16 


16.7586 


-7.020x10" 




—5.555x10" 


2 


-2.158x10" 




-1.733x10" 




8.064x10" 




/ .Uj / X 1U 


2 


2.0 


0.4 


24 


24.9396 


-8.902x10" 




—5.844x10" 


2 


-2.645x10" 




-1.778x10" 




1.122x10" 




y.oou x iu 


2 


2.0 


0.4 


32 


32.9389 


-1.150x10" 




— 6.536x 10" 


2 


-3.267x10" 




-1.896x10" 




1.351x10" 




i . iyj x iu 


1 


1.7 


0.5 








-7.421x10" 


2 


-7.401x10" 


~2" 


-1.815x10" 


1 


-1.810x10" 


1 







u 




1.7 


0.5 


8 


8.4736 


-8.957x10" 


2 


-7.168x 10" 


2 


-2.173x10" 


1 


-1.750x10" 


1 


1.379x10" 


1 


1 \y 1 n - 

l.zjoX IU 


1 


1.7 


0.5 


16 


16.8300 


-1.347x10" 


1 


-6.999x10" 


2 


-3.201x10" 


1 


-1.676x10" 


1 


2.611x10" 


1 


1.5 fo X ID 


1 


1.8 


0.5 








-7.589x10" 


2 


-7.620x10" 


2 


-1.993x10" 


1 


-2.000x10" 


1 







U 




1.8 


0.5 


8 


8.4395 


-8.644x10" 


2 


-6.929x10" 


2 


-2.254x10" 


1 


-1.829x10" 


1 


1.005x10" 


1 


y.u/ox ID 


2 


1.8 


0.5 


16 


16.7776 


-1.175x10" 


1 


-7.210x10" 


2 


-3.004x10" 


1 


-1.880x10" 


1 


1.911x10" 


1 


1. /zox IU 


1 


1.8 


0.5 


24 


24.9413 


-1.678x10" 


1 


-7.395x10" 


2 


-4.158x10" 


1 


-1.881x10" 


1 


2.638x10" 


1 


t toc\/ i n — 
Z. 35DX 10 


1 


1.9 


0.5 








-6.646x10" 


2 


-6.620x10" 


2 


-1.855x10" 


1 


-1.849x10" 


1 







U 




1.9 


0.5 


8 


8.4059 


-7.386x10" 


2 


-6.320x10" 


2 


-2.048x10" 


1 


-1.768x10" 


1 


7.312x10" 


2 


O.j /9X IU 


2 


1.9 


0.5 


16 


16.7233 


-9.572x10" 


2 


—6.551 x 10" 


2 


-2.603x10" 


1 


-1.809x10" 


1 


1.395x10" 


1 


I .ZD j X IU 


1 


1.9 


0.5 


24 


24.8877 


-1.312x10" 


1 


-7.087x10" 


2 


-3.461x10" 


1 


-1.909x10" 


1 


1.937x10" 


1 


1 . / L +H A iu 


1 


1.9 


0.5 


32 


32.8741 


-1.795x10" 


1 


-8.247x10" 


2 


-4.544x10" 


1 


-2.091x10" 


1 


2.320x10" 


1 


2.092x10" 


1 


2.0 


0.5 








-5.987x 10" 


2 


-5.995x10" 


2 


-1.761x10" 




-1.763x 10" 














2.0 


0.5 


8 


8.3750 


-6.516x10" 


2 


-5.918x10" 


2 


-1.906x10" 


1 


-1.738x10" 


1 


5.456x10" 


2 


4.882x10" 


2 


2.0 


0.5 


16 


16.6725 


-8.081x10" 


2 


-5.817x10" 


2 


-2.324x10" 


"1 


-1.694x10" 


1 


1.044x10" 


1 


9.343x10" 


1 


2.0 


0.5 


24 


24.8347 


-1.063x10" 


1 


-6.254x10" 


2 


-2.970x10" 


"1 


-1.776x10" 


1 


1.456x10" 


1 


1.304x10" 


1 


2.0 


0.5 


32 


32.8378 


-1.412x10" 


1 


-6.993x10" 


2 


-3.787x10" 


1 


-1.893x10" 


1 


1.756x10" 


1 


1.576x10" 


1 



Table IV: As in TableslTflandlllll but for different values of eccentricity e; the Teukolsky-based fluxes for E and L z have an accuracy of 10 



15 





a- 

L*inc 




At/M 


A#- 




Al 

V ut "tW 













1.250xl0 6 















5 


5.355510 


1.217xl0 6 


1.949x10" 


1 


4.954x10" 


1 





10 


10.679331 


1.118x 10 s 


3.468x10" 


1 


8.631x10" 


1 





15 


15.943192 


9.574xl0 5 


4.236x10" 




1.019 




o 


20 


21.125167 


7.446 x 10 5 


4.109x 10" 


1 


9.440x 10" 


1 


o 


25 


26.211779 


4.981 xlO 5 


3.158x 10" 


1 


6 860x 10" 

U.Uvv A I V/ 


1 


o 


30 


31 199048 


2.528x 10 5 


1.732x10" 


1 


3.527x 10" 


1 


o 


35 


36.092514 


6.584xl0 4 


4.636x10" 


2 


8.806x10" 


2 


0.1 








1.228x10° 












0.1 


5 


5.351602 


1.198xl0 6 


4.517x10" 


1 


7.766x10" 


1 


0.1 


10 


10.671900 


1.103xl0 6 


6.900x10" 


1 


1.236 




0.1 


15 


15.932962 


9.426 xlO 5 


7.283x10" 


1 


1.344 




0.1 


20 


21.113129 


7.315xl0 5 


6.433x10" 


1 


1.187 




0.1 


25 


26.199088 


4.900 xlO 5 


4.780x10" 


1 


8.547x10" 


1 


0.1 


30 


31.186915 


2.513xl0 5 


2.730x10" 


1 


4.585x10" 


1 


0.1 


35 


36.082095 


6.589xl0 4 


8.385x10" 


2 


1.279x10" 


1 


0.2 








1.173x10° 












0.2 


5 


5.339916 


1.150xl0 6 


1.204 




1.598 




0.2 


10 


10.649670 


1.064 xlO 6 


1.698 




2.331 




0.2 


15 


15.902348 


9.043xl0 5 


1.618 




2.293 




0.2 


20 


21.077081 


6.980 xlO 5 


1.324 




1.900 




0.2 


25 


26.161046 


4.693xl0 5 


9.545x10" 


1 


1.351 




0.2 


30 


31.150481 


2.486 xlO 5 


5.674x10" 


1 


7.711x10" 


1 


0.2 


35 


36.050712 


7.562xl0 4 


2.070x10" 


1 


2.648x10" 


1 


0.3 








1.087 xlO 6 












0.3 


5 


5.320559 


1.069xl0 6 


2.307 




2.788 




0.3 


10 


10.612831 


1.001 xlO 6 


3.256 




4.007 




0.3 


15 


15.851572 


8.454x10 s 


2.984 




3.741 




0.3 


20 


21.017212 


6.483 x10 s 


2.375 




2.998 




0.3 


25 


26.097732 


4.408 x10 s 


1.700 




2.129 




0.3 


30 


31.089639 


2.493 x10 s 


1.040 




1.276 




0.3 


35 


35.997987 


1.108x10 s 


4.626x10" 


1 


5.569x10" 


1 



Table V: Variation in the inclination angles l and #i nc as well as time needed to reach the separatrix for several inspirals through the nearly 
horizon-skimming regime. In all of these cases, the binary's mass ratio was fixed to /i/M = 10~ 6 , the large black hole's spin was fixed to 
a — 0.998A/, and the orbits were begun at p — 1.9M. The time interval At is the total accumulated time it takes for the inspiralling body 
to reach the separatrix (at which time it rapidly plunges into the black hole). The angles A6i nc and At are the total integrated change in these 
inclination angles that we compute. For the e = cases, inspirals are computed using fits to the circular- Teukolsky fluxes of E and L z ; for 
eccentric orbits we use the kludge fluxes d40| l, l |43t and J44t . Notice that A^inc and At are always positive — the inclination angle always 
increases during the inspiral through the nearly horizon-skimming region. The magnitude of this increase never exceeds a few degrees. 



